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We discuss the size of the higher order terms in the NRQCD expansion of the pseudoscalar decay constant. 
Power law divergences in the matrix elements contributing to fps are also investigated. 
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1. Introduction 

When investigating the physics of the heavy 
light mesonic system within the frame work of 
non-relativistic QCD (NRQCD), the action and 
observables are expanded in powers of Aqcd /Mq , 
Mq denoting the heavy quark mass. If the reg- 
ularisation, in our case the lattice cutoff, would 
be removed, divergences will arise in the matrix 
elements and coefficients of the expansion. To 
circumvent this, we keep the cutoff finite and use 
an improved action to minimise residual effects 
on physical results. Large residual divergences 
would become obvious as scaling violations of the 
contributing matrix elements. 

We use 278 quenched 12 3 x 24 configurations 
with j3 = 5.7, generously provided by UKQCD. 
For the light quarks the tadpole improved clover 
action with n = 0.1400 has been used. This is 
slightly lighter than the strange quark mass as 
determined from the K meson, k s \k — 0.1398 [Q. 
At P = 6.0 [0 such a mismatch leads to a shift of 
1% in /ps which is negligible compared to other 
errors in our calculation. From m p one gets a -1 = 
1.103(11)(50) GeV §. For the heavy quarks we 
used a tadpole improved action corrected up to 
order 0(1/M 2 ) in the bare heavy quark mass, 
details on which can be found in 0] . We use heavy 
quark masses in the range 20 > aM > 0.6, which 
ranges from heavy quarks much heavier than the 
b and slightly lighter than the c quark. 

All results presented arc preliminary 
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2. Size of /ps contributions 

The pseudoscalar decay constant /ps is defined 
by the matrix element 



P„/ps := <0|A M |PS> 



(1) 



of a pseudoscalar meson state |PS) and the axial 
vector current A^. In lattice NRQCD A^ has to 
be matched by a set of local operators J 
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M denotes the bare heavy quark mass. The Cj 
are determined in one loop PT |). Since the 
scale q* of the coupling constant a s has not been 
calculated, we will average the final results over 
aq* — 1 and 7T. The difference is treated as a 
systematic error of the pert, expansion. 

In figure [l] we display the size of the matrix ele- 
ment contributing to /ps in different orders of the 
bare heavy quark mass. At the 6-quark the dif- 
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Figure 1. Size of the tree level contribution to 
the pseudoscalar decay constant versus strange 
pseudoscalar meson mass. 



ferent orders in 1 /M are nicely suppressed by one 
order of magnitude with respect to the previous. 
However at the c-quark the O(jrfz) contribution 
has the same size as the 0(4?) part. 

If you look it in detail, this is not that surpris- 
ing. For the term in 0{ i) we observe a suppres- 
sion of 0.36 GeV/M at the B s and 0.29 GeV/M 
at the D s meson with respect to the leading 
term. For the O(^) term the suppression is 
0.26 GeV 2 /M 2 resp. 0.22 GeV 2 /Af 2 . These num- 
bers are compatible with the expectation of a fac- 
tor w Aqcb/M for each order in A, but for the 
O(^p-) term there is apparently a prefactor of ~ 2 
to 3. At the c quark aM < 1, which could have 
some affect on the outcome, requiring further in- 
vestigation. 

3. Power law contributions 

Due to operator mixing, the coefficients c, of 
eqn. (||) would diverge for a — > as O(-^jj), 
P( , % ), ■ ■■J6|, which is unphysical. There- 
fore the divergences have to cancel against simi- 
lar terms arising from the matrix elements of the 
J\\ at - To investigate the size of the unphysical 



Figure 2. O^lat vs strange pseudoscalar meson 
mass. 



part of the matrix elements in comparison to their 
physical, we compare to the results of 0], which 
were obtained at a~ 1 (m p ) = 1.92 GeV. For the 
comparison we study 

0$.t = (OIMJ^JPS), 

= (0|M 2 (jS at + ji 4 j at + jS at )|PS),(4) 

which have a reduced M dependence and a non- 
zero static limit. The result for Ojt\ a t is shown 
in figure 0. No scaling violation is found within 
the quite small error bars. For the ratio of these 
values we obtain 1.06(5) at the B s , which seems 
incompatible with a pure a s /aM power law be- 
haviour. Using a s (aq* = 2) this would lead to a 
ratio of 1.39. Our ratio of the 0^\ at is indepen- 
dent of the quark mass and its error stays below 
0.1 even for the heaviest quark mass. 

For (ft we observe some scaling violation as 
displayed in figure |[ At the B s the result on the 
fine lattice is by a factor of 1.35(7) larger than 
the one from the coarse one. However a power 
law behaviour a s j (aM) 2 would lead to a factor of 
2.4. Again the scaling violations are independent 
of the heavy quark mass. 
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Figure 4. Strange pseudoscalar decay constants 
for two different lattice spacings. 



4. Pseudoscalar decay constant 

In fig. U we show the scaling behaviour of /ps 
itself. The error bars encompass the statistical 
errors and the uncertainty of q* for both graphs, 
since the latter need not to be equal when chang- 
ing the a value. Apart from the heaviest meson 
masses, the error bars on the f3 — 5.7 graph are 
dominated by the uncertainty of q* . The cen- 
tral values differ by 15% ± 10%, which is compat- 
ible with the expected ©((oAqcd) 2 ) discretisa- 
tion correction. Furthermore, a plot of the ratio 
of these results is flat within error bars. 

For the decay constant of the B s we obtain from 
the data set at [3 = 5.7 



f B , = 201(6)(15)(7) MeV. 



(5) 



Here we only give the statistical error, the uncer- 
tainty arising from q* and the determination of 
the bare b quark mass. 

5. Conclusion 

In this talk we discuss the size of the contri- 
butions to fps in a large range of heavy quark 
masses, encompassing the b and the c quark. At 
the b we find higher orders to be nicely sup- 



pressed. Scaling studies find no evidence for large 
power divergent contributions to higher order ma- 
trix elements, which appear to be dominated by 
their physical 0(A QC d/M) resp. C(A| CD /M 2 ) 
term. For the decay constant itself we observe 
reasonable scaling behaviour. 
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